Gauge Invariant Evaluation of Nuclear Polarization with Collective Model 
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The nuclear-polarization (NP) energies with the collective model commonly employed in the NP 
calculations for hydrogenlike heavy ions are found to have serious gauge violations when the ladder 
and cross diagrams only are taken into account. Using the equivalence of charge-current density with 
a schematic microscopic model, the NP energy shifts with the collective model are gauge invariantly 
evaluated for the Isi/a states in ^g|Pb'*^+ and ^||U^^+. 
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High-precision Lamb-shift measurement on high-Z hy- 
drogenUke atoms Q spurred a renewed interest in the 
quantum electrodynamical (QED) calculation of elec- 
tronic atoms. Comparison of theoretical results with ex- 
perimental data allows sensitive tests of QED in strong 
electromagnetic fields ||, |^. In this context, the study 
of the nuclear-polarization (NP) effect becomes impor- 
tant because the NP effect, as a non-QED effect which 
depends on the model used to describe the nuclear dy- 
namics, sets a limit to any high-precision test of QED. 

A relativistic field-theoretical treatment of NP cal- 
culation was presented by Plunien et al. |^ uti- 
lizing the concept of effective photon propagators with 
nuclear-polarization insertions. In these studies, only the 
Coulomb interaction was considered based on the argu- 
ment that the relative magnitude of transverse interac- 
tion is of the order of (f /c)^ and the velocity v associated 
with nuclear dynamics is mainly nonrelativistic. 

The effect of the transverse interaction was studied in 
the Feynman gauge by Yamanaka et al. |7| with the same 
collective model used in ||, ^ for nuclear excitations. 
They found that the transverse contribution is several 
times larger than the Coulomb contribution in heavy 
electronic atoms before the contributions of the positive 
and negative energy states cancel. However, due to the 
nearly complete cancellation between them, the trans- 
verse effects become small and the net effect is destruc- 
tive to the Coulomb contribution in both lsi/2 states of 
As a result, the total NP energy 



20|Pb8i+ and 2||U9i+ 
almost vanishes in ^g|Pb*^+. 
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Recently, the NP effects for hydrogenlike and muonic 
Pb^^+ were calculated in both the Feynman and 
Coulomb gauges, using a microscopic random phase ap- 
proximation (RPA) to describe nuclear excitations ^. 
It was found that, in the hydrogenlike atom, the NP ef- 
fects due to the ladder and cross diagrams have serious 
gauge dependence and inclusion of the seagull diagram 
is indispensable to restore the gauge invariance S. In 



contrast, the magnitude of the seagull collection is a few 
percent effect in the muonic atom, although it improves 
the gauge invariance |^ . 

In the present paper, we report that the nuclear collec- 
tive model employed for hydrogenlike ions in iQ, |[ ^, 0] 
also leads to a large violation of gauge invariance as far 
as the ladder and cross diagrams only are considered. 
Then it is shown, based on the equivalence of the tran- 
sition density of the collective model and a microscopic 
nuclear model with a schematic interaction between nu- 
cleons, that the seagull corrections should also be calcu- 
lated with the collective model in order to obtain gauge 
invariant NP results. The resulting gauge invariant NP 
energy shifts are given for the Isi/2 states in ^glPb*^"*" 
and 2||U91+. 

For spherical nuclei, the Hamiltonian of the small am- 
plitude vibration with multipolarity L is written as 



H, 



^ M Al 



(^LM^LM) 



(1) 



where ttlm are the canonically conjugate momenta to 
the collective coordinates ctLM- The lowest vibrational 
modes are expected to have density variations with no 
radial nodes, which may be referred to as shape oscilla- 
tions. The corresponding charge density operator with 
the multipolarity L is written as 



hit, r) = pL{r) Yl^,jaLMit) 

M 



(2) 



to the lowest order of Q!^^^(t). 

The liquid drop model of Bohr (BM)||l^ is a simple 
model of such a shape oscillation obtained by consider- 
ing deformation of the nuclear radius parameter while 
leaving the surface diffuseness independent of angle: 



LM 



(3) 



where Rq is the nuclear radius parameter of the ground 
state. The radial charge density of BM is given by 
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PUr) 



-Ro-j^go{r), 



(4) 



2 



where go{r) is a charge distribution with spherical sym- 
metry. 

If we assume that under distortion, an element of mass 
moves from Vq to r without alteration of the volume it 
occupies, i.e., the nucleus is composed of an inhomoge- 
neous incompressible fluid, a harmonic vibration of an 
originally spherical surface r = tq in the nucleus is given 

by 



r{n) = ro 



LM 



L-2 



For this model we obtain 



PUr) 



..L-l 



L-2 



dr 



Qoir). 



(5) 



(6) 



This version will be hereafter referred to as the Tassie 
Model (TM) 0. In Eqs.(|) and (|), goir) is usually 
taken to be equal to the ground-state charge distribution. 

In either case, the motion of nuclear matter is assumed 
to be incompressible and irrotational, hence the velocity 
field v{t,r) is given by a velocity potential as v{t,r) ~ 
V<I'(i, r). This implies the nuclear current defined by 
J{r) = gQ{r)v{r) yields the transition multipole density 
of current operator 



M 



Note that the JLL+iif) part does not appear in the tran- 
sition density of current operator given by (|^). 

Therefore, in this kind of collective model, the conti- 
nuity equation of charge gives 



[lAELPLir] 



2L+rdr 



(l )JLL-i{r)] = 0, 



(8) 



where AE'^ is the excitation energy of the surface oscil- 
lation. Hence the transition density of current is given 

by 



JLL-i{r) iAEl 



2L + 1 



x^-^PL{x)dx (9) 



in terms of the transition density of charge. If we assume 
the uniform charge distribution go{r) — goQ{Ro — r), we 
obtain, for both BM and TM, 



oscillation applies to the case of the multipolarity L > 2, 
Eqs. ( p^ and (pT| ) with L = 1 give the transition densities 
of the giant dipole resonance given by the Goldhaber- 
Teller model describing the relative motion of neutrons 
and protons [ p^ . For the monopole vibration, it is also 
possible to construct corresponding charge and current 
densities @, 0]. 

In general, the charge conservation relation between 
the charge and current densities is necessary but not suf- 
ficient for the gauge invariance of the NP calculation. 
Unfortunately, it is practically impossible to construct 
a model that incorporates gauge invariance explicitly in 
terms of the collective variables of the model. However, 
it is possible to evaluate the NP energy gauge invariantly 
with the above collective model as is shown below. 

The NP calculations with the collective model as- 
sume that a single giant resonance with spin multipo- 
larity L saturates the energy- weighted B{EL) strength 
for each isospin. In this respect, let us recall the fact 
that the transition densities of charge to the sum-rule 
saturated levels are given in terms of the ground-state 
charge density |Q. This can be seen as follows. For a 
pair of single-particle operators g{r) = g{r)YLM{^) and 
/(^) = f{f)YLM{^)^ the energy- weighted sum rule can 
be generalized to 



JL{t, r) = JLL-i{r) V YlL_^M^LM{t)- (7) 2J, 



L{L+\) 



2L + 1 
47r 2M 



r^drgo{r)[g'{r)f'{r) 



9{r)f{r)l 
(12) 



where go{r) is the charge distribution of the ground 
state normalized as J r'^drgQ{r) — Z iQ. When a sin- 
gle excited state \ JfMf) saturates the B{EL) strength, 
\JfMf) oc r^YLM\JiMi) , the transition density of charge 
to this state is derived from the sum-rule relation (|l2| ) 
model independently and given by 

If the charge distribution of the ground state is assumed 
to be a uniform distribution with a radius Rq, this be- 
comes 



1 



R, 



L+2 



Sir - Ro), 



(14) 



PLir) = (JfWr^YLU) 



1 



R, 



L+2 



SiRo 



JLL-i{r) - {J}\\r'^YL\\J^) iAEl 



(10) 



2i + 1 r 



L-l 



n2L+ 

Rq 



jOiRo 



(11) 



The transition densities given by (|l^) and (|ri|) have been 
employed in the previous NP calculations for i > 1 ||, 
ItI]. It should be mentioned that, although the surface 



which is equal to the matrix element of the charge density 
operator of the collective model given by (p^). 

On the other hand, it is well known that the schematic 
RPA with a separable interaction 



(15) 



M 



for particle- hole excitations |toz ^) with a degenerate 
particle-hole excitation energy e gives a collective state 
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|LM), which exhausts the energy- weighted sum rule for 
the single particle operator r^Yi^M- 

AEl |(LAf |r^yLA/|0)P - Kmlr^FiMN)!'. (16) 

rni 

where AEl is the excitation energy of \LM) . If the 
ground state is assumed to be a filled major shell of the 
harmonic oscillator potential: 

the particle-hole excitation energy e is taken to be Ihuj 
for 1~ and 2Hiu for 0+ and 2+. The corresponding col- 
lective states exhaust the energy- weighted sum rules, be- 
cause the transition strengths vanish outside these p-h 
excitation spaces. Therefore, the transition densities of 
charge to the collective states of this fictitious nucleus 
are given by ( p^ ) . When the ground-state charge density 
is approximated by a uniform charge density, the transi- 
tion density of charge becomes identical with that of the 
collective model employed in NP calculations for hydro- 
genlike atoms. However, the gauge invariant electromag- 
netic interaction of this schematic microscopic model is 
given by the minimal substitution Pi Pi — SiA to the 
Hamiltonian H — Hfjo+^s- Hence the lowest-order con- 
tributions to NP with this model are given by the three 
Feynman diagrams in Fig. 1, where two photons are ex- 
changed between a bound electron and a nucleus. The 
nuclear vertices are understood to have no diagonal ma- 
trix elements for the ladder and cross diagrams, and no 
nuclear intermediate states for the seagull diagram. It is 
well known that the NP results with this model is gauge 
invariant provided these three diagrams are taken into 
account. Although JLL+i{r) current density appears in 
this model, JLL-iij") dominates in the transition to the 
collective state. 

Thus we can conclude that the gauge invariance of 
the collective model is also guaranteed with the charge- 
current density satisfying the continuity equation (|^), 
provided the contributions from the three diagrams are 
taken into account. It should be noted that the seagull 
contribution is given in terms of the ground-state charge 
distribution and does not depend on the details of the 
model for nuclear excitations. In the actual NP calcula- 
tions 1^, |[ ^, 0] with the collective model, the assumption 
that each nuclear intermediate state saturates the sum- 
rule is not strictly obeyed, because the observed nuclear 
data are used for the low-lying states. However, since 
the gauge violation is serious only in the dipole giant 
resonance, this does not invalidate our arguments as is 
confirmed by the numerical results in the following. 

The formulas to calculate the NP energy shifts due 
to the three diagrams of Fig. 1 were given in for ar- 
bitrary nuclear models. In the present NP calculations 
of the lsi/2 states in hydrogenlike ^gfPb and ^|fU, the 
parameters of the collective model are the same as those 
given in Refs. ^. The same low-lying states and giant 




Electron Nucleus Electron Nucleus Electron Nucleus 
(a) (b) (c) 

FIG. 1: Diagrams contributing to nuclear polarization in low- 
est order; (a) ladder, (b) cross and (c) seagull diagrams. 



TABLE L Nuclear-polarization correction (meV) to the lsi/2 
state of ■^gfPb**^"'' . NP denotes the correction due to the whole 
of the Coulomb and transverse interactions; CNP the correc- 
tion only due to the Coulomb interaction. Energy shifts in 
the parentheses are due to seagull contribution. 



present work Ref. @| Ref. |] 





Feynman(NP) 


Coulomb(NP) 


CNP 


NP 


CNP 


0+ 


-3.3 


(-0.2) 


-3.3 (+0.0) 


-3.3 


-6.6 


-3.3 


1" 


-22.1 


(-42.3) 


-21.5 (-7.3) 


-17.0 


+16.3 


-17.6 


2+ 


-5.8 


(+0.3) 


-5.8 (+0.6) 


-5.8 


-7.0 


-5.8 


3" 


-2.7 


(+0.2) 


-2.8 (+0.2) 


-2.9 


-2.9 


-2.6 


4+ 


-1.0 


(+0.1) 


-1.0 (+0.1) 


-1.1 






5' 


-0.5 


(+0.1) 


-0.6 (+0.0) 


-0.6 






total 


-35.4 


(-41.8) 


-35.0 (-6.4) 


-30.7 


-0.2 


-29.3 



resonances are taken into account. In addition, the con- 
tributions from the 4~ and 5^ giant resonances are also 
calculated in order to see the effects of higher multipoles 
neglected previously. The B{EL) values are adjusted to 
the observed values for low- lying states and the B{EL) 
are estimated through the energy-weighted sum rule for 
giant resonances. 

Tables I and II show the results for the Is 1/2 states in 
^g|Pb^-^+ and ^||U^^+, where the sum of the contribu- 
tions from the three diagrams of Fig.l is given for each 
multipole. The second and the third columns are the re- 
sults including the transverse effects in the Feynman and 
Coulomb gauges, respectively. The values in the paren- 
theses are the contributions from the seagull diagram. 
The NP energy shifts due to the ladder and crossed di- 
agrams only are obtained by subtraction of the seagull 
contributions given in the parentheses. The fourth col- 
umn gives the results of the present Coulomb nuclear 
polarization (CNP). The last two columns are the results 
of the previous calculations. 

The results with the collective model, as with the mi- 
croscopic RPA model also lead to large violations of 
gauge invariance if ladder and crossed diagram contri- 
butions only are considered. The seagull corrections are 
considerable in the 1^ contributions for both of ^glPb*^"*" 
and ^1111^^+. Note that, in the limit of point nucleus, 
which is not unrealistic even for heavy hydrogenlike ions. 
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TABLE II: Nuclear-polarization correction (meV) to the lsi/2 

238 r 
92 ' 



present work 



Ref. 



Ref. § 





Feynman(NP) 


Coulomb (NP) 


CNF 


NP 


CNP 


0+ 


-9.3 (-0.4) 


-9.3 (+0.0) 


-9.3 


-21.5 


-9.5 


1" 


-54.3 (-65.7) 


-52.5 (-3.9) 


-41.6 


-3.8 


-42.4 


2+ 


-131.6 (+0.0) 


-131.7 (+1.6) 


-131.6 


-148.2 


-138.9 


3" 


-6.5 (4-0.3) 


-6.5 (+0.4) 


-6.7 


-7.3 


-6.8 


4+ 


-2.0 (+0.2) 


-2.0 (+0.2) 


-2.1 






5~ 


-1.0 (+0.1) 


-1.0 (+0.1) 


-1.1 






total 


-204.7 (-65.5) 


-203.0 (-1.6) 


-192.4 


-180.8 


-197.6 



the seagull collection occurs only in the dipole mode 
which involves the current density JiQ{r). 

In ^glPb^^+j the contributions from low- lying states 
are about 10% of the total results and the NP energy 
shift is mainly determined by the giant resonance contri- 
butions. The most dominant contribution comes from 
the giant dipole resonance, where a large violation of 
gauge invariance occurs if the seagull contributions in 
the parentheses are neglected: —22 meV becomes +20 
meV and —14 meV in Feynman and Coulomb gauges, 
respectively. The column 5 gives the previous results in 
the Feynman gauge without seagull contributions. The 
differences between the two results in the Feynman gauge 
without seagull contribution come from the accuracy of 
numerical integration over the continuum threshold re- 
gion of electron intermediate states and from the differ- 
ences of the electron wave functions: here we have used 
wave functions in a finite charge distribution, while |^ 
employs point Coulomb solutions. 



In ^g|U^^+, the dominant contribution comes from the 
lowest excited states 2+ with a large B{E2) value. Since 
the transition density of current in the present model 
given by (|ll| ) is proportional to the excitation energy, the 



transverse contribution of the lowest 2+ is negligible due 
to its exceptionally small excitation energy AE2 — 44.9 
keV. Apart from this large Coulomb contribution, the 
contributions from other states show similar tendencies 



as m 



^g2Pb*^+. Namely, the contributions from other 



low-lying states are small compared with the giant reso- 
nance contributions, and a large gauge violation occurs in 
the giant dipole resonance when the seagull contribution 
is omitted. 

To summarize, the transverse effects with the collec- 
tive model are estimated gauge invariantly by inclusion 
of the seagull contribution. The gauge invariance is satis- 
fied to a few percent levels in both ^^^Ph^^+ and ^9211^^+ 
for each of the multipoles separately. Without the seag- 
ull correction, the Feynman gauge in particular does not 
give reliable predictions of NP, although numerical calcu- 
lation in this gauge is easier than in the Coulomb gauge. 
Hence the conclusion of 0| on the transverse effects is no 
longer tenable. The NP energy shifts are — 35.0(— 35.4) 
meV in 20|Pb8i+ and -205(-203) meV in 238u9i+ for 
Coulomb (Feynman) gauge. The net transverse effect is 
about 14 ^ 15% of the Coulomb energy shift of —30.7 
meV in ^g|Pb^"'^+. This is similar to the conclusion of 
the microscopic model [||, and should be compared with 
the transverse effect of the lsi/2 state in muonic ^gfPb, 
which is about 6 % of the Coulomb contribution § . The 
agreement between the two models provides stability of 
the prediction of the NP effects with respect to the choice 
of the nuclear models. The percentage of the transverse 
effect in the total shift in ^gfU^^^ is reduced to about 
6% of the Coulomb effect due to the dominant Coulomb 
contribution from the lowest 2+ state. 
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generous support and useful discussions during our re- 
search on NP effects. They appreciate Drs. N. Yamanaka 
and A. Ichimura for collaboration on the NP effects with 
the collective model, which motivated the present work. 
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